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ABSTRACT 

High-resolution hydrodynamics simulations of the Lya forest in cold dark matter 
dominated cosmologies appear to predict line-widths that are substantially narrower 
than those observed. Here we point out that Compton heating of the intergalactic gas 
by the hard X-ray background (XRB), an effect neglected in all previous investigations, 
may help to resolve this discrepancy. The rate of gain in thermal energy by Compton 

Cn| ' scattering will dominate over the energy input from hydrogen photoionization if the 

XRB energy density is ~ 0.2 x/ (e) times higher than the energy density of the UV 

Q^ ■ background at a given epoch, where x is the hydrogen neutral fraction in units of 10 -6 

and (e) is the mean X-ray photon energy in units of m e c 2 . The numerical integration of 
the time-dependent rate equations shows that the intergalactic medium approaches a 
temperature of about 1.5 x 10 4 K at z > 3 in popular models for the redshift evolution 
of the extragalactic background radiation. The importance of Compton heating can be 
tested experimentally by measuring the Lya line-width distribution as a function of 
redshift, thus the Lya forest may provide a useful probe of the evolution of the XRB 

S^ ■ at high redshifts. 

Subject headings: cosmology: theory - diffuse radiation - intergalactic medium - quasars: 
absorption lines - X-rays: general 



1. Introduction 

Neutral hydrogen in the intergalactic medium (IGM) produces a plethora of Lya absorption 
lines in the spectra of high redshift quasars. Numerical N-body/hydrodynamics simulations of 
structure formation in the IGM, within the framework of cold dark matter (CDM) dominated cos- 
mologies (Cen et al. 1994; Zhang, Anninos, & Norman 1995; Miralda-Escude et al. 1996; Hernquist 
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et al. 1996; Zhang et al. 1998; Theuns et al. 1998a, b; Dave et al. 1999; Bryan et al. 1999), have 
recently provided a definite picture for the origin of the Lya forest, one of an interconnected net- 
work of sheets and filaments with virialized systems (halos) located at their points of intersection. 
At z > 2, the lowest column density absorbers (Nm ~ 10 12 cm~ 2 ) arise in this context in the un- 
derdense (pb/pb < 1) minivoids between the filaments, absorbers that give rise to moderate column 
density lines (iVm ^ 10 14 cm -2 ) correspond to modestly overdense (1 < Pb/pb < 5) filaments, and 
high column density absorbers (iVni ^ 10 15 cm -2 ) arise in highly overdense (pb/pb > 10) structures 
at the intersection of filaments (Zhang et al. 1998). Most simulations to date assume that the 
IGM is photoionized and photoheated by a UV background close to that inferred from quasars (as 
computed by Haardt & Madau 1996, hereafter HM96). The overall normalization of the column 
density distribution depends approximately on the parameter i7 2 /i 3 /T, where T is the hydrogen 
ionization rate, h is the Hubble constant in units of 100 kms -1 Mpc -1 , and f^ is the baryon den- 
sity parameter. The width of the lines, as measured by the b parameter of a Voigt profile, is set by 
thermal broadening, peculiar velocities, and Hubble expansion across the filaments. 

While the first simulations showed good agreement with the observed line statistics, recent 
detailed studies at higher numerical resolution have revealed a serious conflict with the data: the 
models predict median values for the b distribution of about 20 kins -1 at z = 3, compared with the 
observed median of 30 kms -1 (Theuns et al. 1998b; Bryan et al. 1999). One possible way to restore 
the agreement with observations is to increase the temperature of the IGM, adding to the thermal 
broadening. As Theuns et al. (1999) have pointed out, this may be achieved by increasing the baryon 
density w& = Qbh 2 , or by heating the gas for longer by increasing the age of the universe (e.g. by 
lowering the total matter density 0, m ). However, the gas temperature at late times is determined 
(approximately) by the balance between adiabatic cooling and photoelectric heating and varies as 
T oc (f^/i/y^m) 1 ' 1 ' 7 (Hui & Gnedin 1997). Thus the needed increase in temperature (a factor of 
2) is difficult to accomplish without conflicting with big bang nucleosynthesis for reasonable values 
of £l m and h (Theuns et al. 1999 and § 3 of this paper). Delaying the reionization of helium may 
also help (Haehnelt & Steinmetz 1998), but this appears to boost the temperature by only a small 
factor. Other mechanisms that could increase the line-width include photoelectric dust heating 
(Nath, Sethi, & Shchekinov 1999) and radiative transfer effects during the reionization of helium 
(Miralda-Escude & Rees 1994; Abel k Haehnelt 1999). 

In this Letter we point out that a potentially important heating term has been neglected in 
all previous numerical studies, namely Compton heating of electrons by hard X-ray background 
photons. We show that this effect will dominate over photoelectric heating of hydrogen at redshifts 
> 2 in popular models for the redshift evolution of the extragalactic background radiation and that 
it will drive the IGM towards a temperature of about 1.5 x 10 4 K at z > 3. 



2. Compton heating of the IGM 

In an IGM that interacts solely with a background of UV ionizing photons (assumed to originate 
from quasars and/or star-forming galaxies), the dominant heat input is due to photoelectric heating 
of hydrogen and helium atoms. The energy input from hydrogen photoionization is given by 

/■OO j 

#H = 47rnHi/ dv-r—an(hpv-hpVL), (1) 

Ju L hpv 

where nni is the mean neutral hydrogen density, J v oc (u jvL)~ a the specific intensity of the UV 
radiation background, cth = <Jl(v/vl)~ 3 the photoionization cross-section, and hpu^ the energy of 
the Lyman edge. Equation (|l|) can be integrated to give 

Hu = n m a L c U vv t — -^ — — r (2) 

(a + 2)(a + 3) 

where U\jy is the (proper) energy density of the metagalactic flux. In a quasar-dominated back- 
ground, photoheating by He n will be larger than Hr by a factor 1.5-2. The main cooling processes 
are Compton cooling against the microwave background at high redshifts and adiabatic cooling at 
redshifts < 2. 

Allow now the intergalactic gas to exchange energy with the XRB through Compton heating 

H c = 4vr n e — ^r / dv J v (h P v - 4kT), (3) 

m e c z Jo 

where n e is the mean electron density, J v the specific intensity of the XRB, and ax the Thomson 

cross-section. Energy is transferred from photons to electrons when hpv 3> 4kT, and in this limit 

equation (||) becomes 

H c = a T n e c { -^4u x , (4) 

m e c z 

where {hpv) is the mean photon energy and Ux the energy density of the XRB. For photons above 

100 keV relativistic effects become significant, and Compton scattering takes place in the Klein- 

Nishina regime. The relativistic generalization of equation (||) can be written as (cf. Blumenthal 

1974) 
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(for e 3> 4kT/m e c 2 ), where e is the photon energy in units of m e c 2 . All the energy of the primary 
electrons will be deposited into the medium as heat (through Coulomb collisions), since energy 
losses due to collisional excitation and ionization of H and He are negligible in such a highly ionized 
IGM (e.g. Shull & Van Steenberg 1985). 

The XRB spectrum is characterized by a narrow peak in spectral power, EJg, at 30 keV; a 
simple analytical fit to the observed energy flux in units of keVcm~ 2 s _1 keV -1 sr _1 yields 

-0.29 „,,„ ( -^kcV 



Je = 7.7E-^ exp ( -^ ) (6) 



from 3 to 60 keV (Boldt 1987). At higher energies, the HEAO 1 A-4 data can be fitted with a 
power-law, 

/ -p \ -1.75 

* = °- 26 (w) m 

from 80 to 400 keV (Kinzer et al. 1997). The preliminary diffuse 7-ray spectrum from 800 keV to 
30 MeV measured with COMPTEL (Kappadath et al. 1996) obeys another power-law, 

J E = 0.15 (^ kc v/100)- L38 keVcm" 2 s" 1 keV" 1 sr" 1 (8) 

(Kribs & Rothstein 1997). The total energy density in the 0.003-30 MeV band is U x = 6.3 x 
10" 5 eVcm" 3 . The mean, "effective" (as it includes a large correction due to Klein-Nishina effects) 
photon energy is (hpf) = 31.8 keV. Note that, because of the decline at high energies of the Klein- 
Nishina cross-section, the (X+7)-ray heating rate is numerically equivalent to the rate one would 
derive by assuming pure Thomson scattering of photons in the range 3-180 keV only. With the 
adopted spectrum, about half of the total energy exchange rate at the present-epoch is provided 
by photons above 110 keV and about 14% is contributed by photons above 1 MeV. 

We can now estimate the relative importance of Compton versus photoelectric heating of 
hydrogen. The QSO contribution to the local metagalactic flux at the Lyman edge is Jl ~ 10 -23 
ergs cm" 2 s _1 Hz" 1 sr" 1 (HM96). Assuming a spectral index of a = 1.3, from equations (|l]) and (|3|) 
we get 

— - « 6.8 x lO" 7 ^- (9) 

H H n H i 

at z = 0. The equation of ionization equilibrium for an IGM with mean hydrogen density 
nn = 1-6 x 10" 7 (u^/0.019) cm" 3 and temperature T = 2 x 10 3 K gives n e /nHi ~ 1.5 x 10 5 , 
hence Hq ~ 0.1 i^n- At the present-epoch, photoelectric heating is then the dominant contribu- 
tion. At z = 2, the UV background from optically-selected quasars increases to Jl ~ 3 x 10" 22 
ergs cm" 2 s _1 Hz" 1 sr" 1 (HM96), the gas temperature reaches 10 K or higher, and ionization equi- 
librium now yields n e /nni ~ 5 x 10 5 . If the bulk of the XRB is produced at high redshifts, Ux 
increases as (1 + z) 4 , and the rate of gain in thermal energy per baryon due to Compton scattering 
goes up as (1 + z) 4+,a , where /3 ~ 1/3 is less than unity because Klein-Nishina corrections become 
more important at high redshifts. The same exercise then yields He ~ 1.3 Hn at z = 2. 

The relative strength of the two competing heating terms will obviously depend on the redshift 
evolution of ionizing and X-ray extragalactic radiation. Two cosmological effects may tend to make 
Compton scattering of hard X-ray photons the dominant source of heating at high redshifts in a 
highly ionized, low density IGM: 

• at fixed ionization state and energy density ratio Ux/U\jy, the ratio Hc/Hr increases with 
(hpu) as (1 + z) 13 , where /3 w 1/3 for z < 3; 

• at fixed ratio between the quasar UV and X-ray volume emissivities, the radiation energy 
density Ux increases with redshift relative to U\jy, as the mean free path of photons at 



912 A becomes very small, sources at higher redshifts are severely attenuated, and the UV 
background is largely "local" beyond a redshift of z > 2 (e.g. Madau 1992). 

In the absence of any other cooling or heating mechanisms, Compton cooling off cosmic microwave 
background (CMB) photons and Compton heating by the XRB will drive the thermal state of the 
IGM towards the Compton temperature 

_ (h P v) U x , n v 

independently of the electron density. At z = 0, Tq ~ 2.2 x 10 K. 



3. Non equilibrium models 

In this section we will calculate numerically the thermal history of a uniform IGM of primordial 
composition within a cosmological context. The code that we use includes the relevant cooling and 
heating processes (as described in Theuns et al. 1998b, and references therein) and follows the 
non-equilibrium evolution of hydrogen and helium ionic species. The gas is allowed to interact 
with the CMB through Compton cooling, with the XRB through Compton heating, and with the 
time-dependent QSO ionizing background as computed by HM96. 

It is widely believed that the diffuse XRB is due to the summed emission from AGNs (see 
Fabian &; Barcons 1992 for a review). There is insufficient information, however, to uniquely infer 
from the data the evolution with cosmic time of the hard X-ray volume emissivity. According to 
the findings of Boyle et al. (1994) from ROSAT data (see also Page et al. 1996), the model that 
best fits the cosmological evolution of soft X-ray selected AGNs in an Einstein-deSitter universe 
is one in which the comoving volume emissivity increases as (1 + z) 315 for z < 1.6 (in qualitative 
agreement with the evolution of the QSO optical luminosity function), and remains roughly constant 
at higher redshifts, up to z ~ 4 (Miyaji, Hasinger, & Schmidt 1998). Assuming that these simple 
prescriptions also describe the evolution of hard X-ray sources, the energy density of the XRB then 
increases approximately as Ux ex (1 + z) 4 at low redshifts, and as (1 + z) 2 ' 5 at z > 1.6. To account 
for the existence of a high redshift cut-off, we will parameterize the evolution of the X-ray energy 
density as 

U x (z) = U X (0){1 + zf exp(-z 2 /^ 2 )- (11) 

Figure 1 shows the ensuing Compton heating rate per baryon as a function of redshift for z c = 5, 
compared with the photoionization heating rate from optically-selected quasars. Figure 2 shows 
the thermal history of an IGM with w& = 0.019 (Buries and Tytler 1998) for two cosmological 
models: an Einstein-deSitter universe and a spatially flat universe dominated by a cosmological 
constant with parameters close to those inferred from Type la supernovae and cosmic microwave 
background anisotropies (see Efstathiou et al. 1999, and references therein). In both cases, we 
assume a photoionizing background as computed by HM96 and a Hubble constant of h = 0.65 (see 



Freedman et al. 1998). The dotted lines show the temperature evolution if the X-ray background 
is ignored, as has been assumed in all previous work. The solid lines include Compton heating by 
the XRB according to the model of equation (|ll|) . Including the X-ray background with a redshift 
cut-off parameter of z c = 5 causes the temperature of the IGM to rise to 1.4 x 10 4 K at z = 4 and 
1.1 x 10 4 K at z = 3. The IGM temperatures are even higher in the A-dominated model. These 
results are extremely insensitive to the amplitude of the photoionizing background and relatively 
weakly dependent on u>& and h. They are, however, sensitive to the assumed redshift evolution of 
the X-ray background and to the hard X-ray spectrum at photon energies above 50 keV, which is 
poorly known but dominate the heating rate. 

The results of Figure 2 and the numerical simulations of Theuns et al. (1999) suggest that 
heating of the IGM by the hard XRB may well be able to account for the 6-parameter distribution 
of the Lya forest. Furthermore, our results suggest that the observed 6-parameter distribution 
may provide a useful probe of the evolution of X-ray diffuse radiation at high redshift. Compton 
heating will also tend to flatten the equation of state of the intergalactic gas, since the Compton 
heating per unit mass is independent of density. However, for our model with z c = 5, this effect is 
small, even for low density regions. 



4. Discussion 

We have shown that Compton heating of electrons by hard X-ray background photons provides 
an important heating source of the IGM. For plausible parameters, heating by the XRB will be 
comparable to photoelectric heating at redshifts ^ 2. This effect has been ignored in most numerical 
simulations of the Lya forest, yet it can have a significant impact on the temperature of the IGM, 
raising it to temperatures as high as ~ 1.5 x 10 4 K at redshifts z > 3. 

The increase in IGM temperature caused by this effect may explain why the Lya 6-parameter 
distributions derived from numerical simulations (Theuns et al. 1998b; Bryan et al. 1999) are 
narrower than those observed. Theuns et al. (1999) have previously suggested that the temperature 
of the IGM can be raised by invoking a high baryon density in a universe with low total matter 
density. However, the density required, w& > 0.025, is well above the value implied by primordial 
nucleosynthesis and the measurements of the deuterium abundances from quasar spectra (cob = 
0.0193 ± 0.0014, Buries and Tytler 1998). Compton heating by hard X-ray photons may help to 
resolve this discrepancy. The origin and main production epoch of the XRB are still largely unsolved 
problems. The ^-parameter distribution of the Lya lines could provide important information on 
the evolution of X-ray background radiation at high redshifts. 

Finally, it is interesting to note that, as another consequence of Compton heating of the 
IGM, there may exist an X-ray proximity effect, in analogy with the usual UV proximity effect 
(e.g. Bechtold 1994 and references therein). This could result in a hotter IGM, and consequently 
broader Lya lines, in the immediate vicinity of a long-lived X-ray emitting quasar (lifetime > the 



Compton heating timescale). 
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Fig. 1. — Compton (dashed line) and photoionization (solid line: H only; dash-dotted line: H+ He) heating 
rates per baryon as a function of redshift. The rates have been derived assuming ionization equilibrium in an 
IGM with ojf, = 0.019, the HM96 UV radiation field, the XRB evolution described in the text (with z c = 5), 
and the corresponding temperature evolution shown in Figure 2b. 
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Fig. 2. — Thermal history of a uniform intergalactic medium for two spatially flat cosmologies. The dotted 
lines show the temperature evolution when the X-ray background is ignored and the only heating source is 
the quasar UV background as computed by HM96. The solid lines show the temperature evolution when the 
XRB is included with redshift evolution as given by equation (|l|) for z c = 5. We have assumed u>b = 0.019 
and h — 0.65. 



